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Abstract. We provide exact and asymptotic formulae for the number of unrestricted, 
respectively indecomposable, d-dimensional matrices where the sum of all matrix en- 
tries with one coordinate fixed equals 2. 

^ ! 1. Introduction 

^ ; We begin by recalling the notion of a magic matria^ this is a square matrix m = 

{ m i.j)i<i.j<n with non-negative integral entries such that all row and column sums are 
equal to the same non-negative integer. If this non-negative integer is s, then we call 
such a matrix s-magic. The enumeration of s-magic squares has a long history, going 
back at least to MacMahon [151 §404-419]. A good account of the enumerative theory 
of magic squares can be found in [181 Sec. 4.6], with many pointers to further literature. 
For more recent work, see for instance [H [8] . 

Let [n] denote the standard n-set {1,2, ... ,n}. There are two obvious ways of gen- 
eralising s-magic matrices to higher dimensions: 

(Gl) All line sums are equal. Given a positive integer d, a <i-dimensional matrix 
m : [n] d — > No (where No denotes the set of non- negative integers) is called 
s-magic if 

^2 rn(u} 1 ,u} 2 ,...,u)d) = s (1.1) 

0Ji&[n] 

for all fixed . . . , Wi-i, Wj+i, . ..,^6 [n], and alii = 1, 2, . . . , d. 



2000 Mathematics Subject Classification. Primary 05A15; Secondary 05A16 05A19 05B15. 

Key words and phrases, higher-dimensional magic matrices, labelled combinatorial structures, mul- 
tisort species, exponential principle. 

■'■Research partially supported by the Austrian Science Foundation FWF, grants Z130-N13 and 
S9607-N13, the latter in the framework of the National Research Network "Analytic Combinatorics 
and Probabilistic Number Theory" . 

1 Strictly speaking, the correct term here would be "s-semi-magic," since we do not require diagonals 
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we prefer the term "magic" for the sake of brevity. 
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(G2) All hyperplane sums are equal. Given a positive integer d, a d- dimensional 
matrix m : [n] d — > Nq is called s-magic if 

m(wi, w 2 , • • • , w<f) = s (1.2) 

wi,...,Wi_i,w< + i,...,a;d6[n] 

for all fixed Uj € [n], and all i = 1, 2, . . . , d. 

Generalisation (Gl) appears already in the literature, see e.g. [HH]. For d = 3 and 
s = 1, these objects are equivalent to Latin squares counted up to isotopy: the roles 
of rows, columns, and symbols of the corresponding Latin square are played by the 
first, second, and third coordinate, respectively, and the entry in position (wi,u; 2 ) of 
the Latin square is U3 if and only if m(ui, 0J2, 0J3) = 1. 

Generalisation (G2) appears in the literature (in more general form) as contingency 
tables in statistics; there are Markov chain methods for approximate counting of these, 
as well as some remarkable asymptotic estimates, see [TTJ [HI EH EE HU]. Indeed, these 
results suggest that the counting problem for (G2) is much easier than for (Gl). (We 
are grateful to a referee for this information and the references.) 

The present note focusses on the second generalisation. Hence, from now on, when- 
ever we use the term "s-magic," this is understood in the sense of (G2). 

Counting higher- dimensional magic matrices is made more difficult (than the already 
difficult case of 2-dimensional magic matrices) by the fact that the analogue of Birkhoff 's 
Theorem (cf. [5] or [2[ Corollary 8.40]; it says that any 2-dimensional s-magic matrix 
can be decomposed in a sum of permutation matrices, that is, 1-magic matrices) fails for 
them. For example, the 3-dimensional 2-magic matrix with ones in positions (1, 1, 1), 
(1, 2, 3), (2, 1, 2), (2, 2, 1), (3, 3, 2) and (3, 3, 3) is not the sum of two 1-magic matrices. 

As we demonstrate in this note, it is however possible to count the 2-magic matrices 
of any dimension. Our first result is a recurrence relation for the number u n (d) of 
indecomposable ^-dimensional 2-magic matrices of size n (see Corollary [3] in Section H]) . 
This recurrence is used in Proposition H] to derive, for fixed d > 3, an asymptotic 
formula for u n (d). In order to go from indecomposable matrices to unrestricted ones, 
we observe that the <i-dimensional 2-magic matrices may be viewed as a <i-sort species in 
the sense of Joyal [2] which obeys the (d-sort) exponential principle. Let w n (d) denote 
the number of all d- dimensional 2-magic matrices of size n. The exponential principle 
can then be applied to relate the numbers w n (d) to the numbers u n (d), see (13.51) (for 
d = 2) and (16. ip (for d > 2). This relation is used in Theorem [5] to find, for fixed d > 3, 
an asymptotic estimate for the numbers w n (d) as well. Exact and asymptotic formulae 
for u n (d) and w n (d) for d — 2 are presented in Section [3j We remark in passing that 
a simple counting argument shows that the obvious interpretation of the matrices in 
Generalisation (Gl) as a d-sort species does not satisfy the exponential principle, not 
even under the — in a sense — minimal axiomatics of [7]. 

2. Indecomposable 2-magic matrices and fixed-point-free involutions 

A d- dimensional matrix m : [n] d — > No is called decomposable, if there exist non-empty 
subsets BP, b£\b[ 2 \ Bf\ B[ d) ,B { 2 d) of [n] with 

B? II B« = Bf ] II Bf = ■■■ = B[ d) II B {d) = [n] 

(II denoting disjoint union) and 

|£«| = |£f)| = ...= | J B( d )| ) 
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such that m(u)i,u)2, • • • , uid) 7^ only if either 

(wi, u 2 , ■ ■ ■ , Wd) e b{ 1} x ) x • • • x B 



or 




of the lines of the matrix such that m attains a block form.) The integer n is called the 

size of m. 

Let u n (d) denote the number of indecomposable <i-dimensional 2- magic matrices of 
size n. Note that an indecomposable 2- magic matrix with an entry 2 has size 1. So it 
is enough to consider zero-one matrices. 

The purpose of this section is to relate the numbers u n (d) to another sequence of 
numbers v n (d) counting certain tuples of fixed-point-free involutions on a set with 2n 
elements. More precisely, let 



be the standard fixed-point-free involution on the set [2n]. Then we define v n (d) to be 
the number of choices of d— 1 fixed-point-free involutions t 2 , ■ ■ ■ , t^ on [2n] such that the 
group G = (ti, t 2 , ■ ■ ■ , t d ) generated by t 1 , t 2 , . . . , t d is transitive. (For example, when 
n = 2, there are just three fixed-point-free involutions on {1,2,3,4}, viz., (1,2)(3,4), 
(1,3)(2,4) and (1,4)(2,3), any two of which generate a transitive group. So v 2 (d) = 



Proof. Let m be an indecomposable <i-dimensional 2-magic matrix of size n, where 
n > 1. Then m is a zero-one matrix, and it contains 2n entries equal to 1, the rest 
being zero. Number the positions of the l's in m from 1 to 2n in such a way that the 
positions with first coordinate j are numbers 2j — 1 and 2j for j = 1, . . . , n. (There are 
2" ways to do this, since for each j we can choose arbitrarily which of the two l's has 
number j — 1.) Then, for i = 1, . . . , d, let tj be the fixed-point-free involution whose 
cycles are the pairs of numbers in {1, . . . , 2n} indexing positions of l's with the same 
i-th coordinate. Note that t\ is the involution defined in (12. ip . 

We claim that the subgroup G of S 2n generated by ti, . . . , td is transitive if and only 
if the matrix m is indecomposable. For this, note that the l's whose labels belong to 
a cycle of ti have the same z-th coordinate. So, if m is decomposable, and the 1 with 
label 1 belongs to x • • ■ x B[ d \ then an easy induction shows that any 1 whose 
label is in the same orbit belongs to this set, so that G is intransitive. Conversely, if 
G is intransitive, then the coordinates of the l's whose labels belong to a G-orbit give 
rise to a decomposition of m. 

So each matrix gives rise to 2 n such (i-tuples of involutions. Thus, the number of 
pairs consisting of a matrix and a corresponding sequence of permutations is 2 n u n (d). 

For instance, the example of a matrix failing the analogue of Birkhoff's Theorem 
given in the Introduction, with the entries numbered in the order given, produces the 
three permutations (1, 2)(3, 4)(5, 6), (1, 3)(2, 4)(5, 6) and (1, 4)(2, 6)(3, 5). 

2 We warn the reader that for d — 2 this does not reduce to the notion of decomposability of matrices 
in linear algebra since there rows and columns are reordered by the same permutation. Yet another 
definition of indecomposability occurs in [I] . 



h = (l,2)(3,4)---(2n-l,2n) 



(2.1) 




v n (d). 



(2.2) 
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Conversely, let t±, . . . , td be fixed-point-free involutions on the set {1, . . . , 2n} which 
generate a transitive group, where t\ is the standard involution defined in (12. ip . Number 
the cycles of each ti from 1 to n such that the cycle (2j — 1, 2j) of ti has number j. (There 
are (n!)^" 1 such numberings.) Now construct a d- dimensional matrix m as follows: for 
k = 1, . . . , 2n, if k lies in cycle number u>i of ti, then m(ui,U2, ■ ■ ■ , ojd) = 1; all other 
entries are zero. Then m is 2-magic. Consequently, each sequence of permutations 
gives rise to (n!) d_1 matrices; and the number of pairs consisting of a matrix and a 
corresponding sequence of permutations equals (n!) d_1 v n (d). 

Comparing these two expressions, we obtain ( I2.2p . as required. □ 

Remark. We note that u\(d) = v\(d) = 1 for all d. Hence, Formula ( 12. 2 p is false for 
n = 1. 



3. Computation of u n {2) and w n (2) 

The number w n (2) of 2-dimensional 2-magic matrices of size n has been addressed 
earlier by Anand, Dumir and Gupta in [21 Sec. 8.1]. They found the generating function 
formula 

E-^T^Hl-^V/ 2 . (3.1) 



n>0 



This gives the explicit formula 



, x A f2k\ (n\f . N 



Singularity analysis (cf. [121 Ch. VI]) applied to ( 13. ip then yields the asymptotic formula 

~ ( „ / 1 



w n (2) = (n!) 2 y — ^1 + C l-j J , as n ->• oo. (3.3) 

The number w n (2) of indecomposable 2-dimensional 2-magic matrices of size n can 
also be computed explicitly. One way is to observe that, by Birkhoff's Theorem (cf. 
[5] or [21 Corollary 8.40]), a 2-magic matrix m is the sum of two permutation matrices, 
say pi and p 2 - If m is indecomposable, then the pair {pi,P2} is uniquely determined. 
Premultiplying by Pi , we obtain a situation where p% is the identity; indecomposability 
forces P2 to be the permutation matrix corresponding to a cyclic permutation, since a 
cycle of P2 not containing all points would provide a decomposition of m. So there are 
n\ (n — 1)! choices for (pi,P2), and half this many choices for m (assuming, as we may, 
that n > 1). Note that this formula gives half the correct number for n — 1. So we 
have 

un{2) = {\ , N * n = 1 ' (3.4) 
V ; |fn!(n-l)!, ifn>l. V ; 

Alternatively, we may observe that 2-dimensional 2-magic matrices may be seen as a 
2-sort species in the sense of Joyal [Hj (see also [HI Def. 4 on p. 102]), with the row 
indices and the column indices forming the two sets on which the functor defining the 
species operates. Hence, by the exponential principle for 2-sort species [TH Prop. 20] 
(see also [HI Sec. 2.4]), we have 

J2 w ^j^y = ex p J2 u ^j^ ■ ( 3 - 5 ) 

n>0 V \n>l V '' J 
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Combining this with (13. ip . we find that 

I>"< 2 >PF = f4 log (rb)- 

n>l K 1 v 7 

Extraction of the coefficient of z n then leads (again) to (13. 4p . 

4. A RECURRENCE RELATION FOR V n (d) 

In this section we prove a recurrence relation for the numbers v n (d) (see Section [2] for 
their definition). By Lemma [U this affords as well a recurrence relation for the numbers 
u n (d). 

Proposition 2. The numbers v n (d) satisfy v\(d) = 1 and 

E (l Z !) ^ 2n ~ 2k - Vk ^ = ^ 2n - n > 1. (4.1) 

fc=i ^ ' 

Here, (2n — 1)!! = 1 ■ 3 • 5 • ■ ■ (2n — 1) is the product of the first n odd positive integers 
forn > 0, and, by convention, (— 1)!! = 1. 

Proof. Recall that (2n — 1)!! is the number of fixed-point-free involutions on a set of 
size 2n. (This is a special case of the general formula 



nl 



for the number of permutations in S n with a« cycles of length i for i = 1, . . . , n.) The 
number of choices of involutions ti,t 2 , ■ ■ ■ , td, where t 1 is as in (12. ip . such that the orbit 
containing 1 of the group they generate has size 2k is 



n — 1 
fc-1 



((2n-2fc- l)!!)^ 1 ^^) 



since we can choose in order 

(i) k — 1 of the n — l cycles of t\ other than (1,2) such that the elements not fixed 
by all of these k — 1 transpositions together with {1,2} form the desired orbit, 
O say; 

(ii) d—1 fixed-point-free involutions on O which, together with the restriction of t\ 
to O, generate a transitive group; 

(iii) d—1 arbitrary fixed-point-free involutions on the complement of O. 

Summing these values shows that the numbers v n (d) satisfy the desired recurrence. □ 
Corollary 3. For all integers d > 1, the numbers u n (d) satisfy Ui(d) = 1 and 

(( 2„ - 3)!o« + ± g : ;) f h -;-" 11 )" = «», - did", 



n > 1. 



Remarks. (1) In the case d = 2, we have seen in (13.41) that u n (2) = n\(n — l)!/2 for 
n > 1, so that 

Vn (2) = 2"- 1 (n-l)! = (2n-2)!!, 
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d 




n = 1 


n = 2 


n = 3 


n = 4 


n = 5 


n = 6 


2 


indec 


1 


1 


6 


72 


1440 


43200 




0-1 





1 


6 


90 


2040 


67950 




all 


1 


3 


21 


282 


6210 


202410 


3 


indec 


1 


8 


900 


359424 


370828800 


820150272000 




0-1 





8 


900 


366336 


378028800 


833156928000 




all 


1 


12 


1152 


431424 


427723200 


920031955200 



Table 1. Indecomposable, zero-one and arbitrary d- dimensional 2-magic 
matrices of size n 



where (2n — 2)!! is the product of the even integers up to 2n — 2 (with 0!! = 1 by 
convention). Substituting this in (I4.ip . we have proved the somewhat curious looking 
identity 

^ (k ~ l) ^ ~ 2k ~ ( 2k ~ 2 ) !! = ( 2n ~ 
k=i ^ ' 

for n > 1. 

We remark that this identity has an interpretation in terms of hypergeometric func- 
tions, for which we refer to [16], in particular, (1.7.7), Appendix (III. 4). The left-hand 
side is 

2"" 1 (1/2)^ ■ a F a "^ n + + V; 
and the identity is an instance of the Chu-Vandermonde identity. 



(2) For d > 2, we have not been able to solve the recurrence explicitly. However, it 
is easy to calculate terms in the sequences, and we can describe their asymptotics (see 
Sections and [6]) • 

Table [1] gives counts of all indecomposable matrices, all zero-one matrices, and all 
non-negative integer matrices, with dimension d and hyperplane sums 2. The sequences 
for d = 2 are numbers A010796, A001499, and A000681 in the On-Line Encyclopedia 
of Integer Sequences [H]. For d = 3, they are A112578, A112579 and A112580. 

5. Asymptotics of the numbers u n (d) 

This section provides the preparation for the determination of the asymptotics of the 
numbers w n (d) for d > 3 in the next section. Our goal here is to establish an asymptotic 
estimate for the sequence u n (d) with fixed d > 3. 

Proposition 4. For fixed d > 3, we have 

u n {d) ~ 2- dn ((2n)\) d -\ asn^oo. 

Proof. By Lemma (TJ we have u n {d) = (n\) d ~ l v n (d) /2 n for n > 1, so it suffices to show 
that 

v n {d)~{{2n-l)\\) d -\ 

We will use the estimates 

. 2 n n 1 
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for n > 1. With c n = 2 n n\/{2n - 1)!!, we have c n+l /c n = {2n + 2)/{2n + 1), and 
both inequalities are easily proved by induction. From these estimates, we obtain the 
inequality 

(2n - 1)!! > fn\ Hk + l)(n - k + 1) \ 1/2 ^ (g ^ 



(2Jfe - 1)!! (2n - 2Jfe - 1)!! - VJv \ n 
To simplify our formulae, we denote the left-hand side of this inequality by 



Now, by Proposition 14. 1\ v n (d) satisfies the recurrence 



n_1 / — 1 \ 

v n (d) = ((2n - l)!!)'- 1 " E U - 1 ((2n ~ 2fc ~ n > L 

fc=i ^ ' 

Clearly u„(d) < ((2n - l)!!)^ 1 . We show that v n (d) > {(2n - l)!!) 4 *" 1 ^ - 0(l/n)), an 
estimate which, in view of the above recurrence, follows if we can show that 

-=£(;:!) (0)-"-e 

Using (15. ip . we have 

n-2 / i\ / \ -W-l) / \ (d-l)/2 



n s --^ fn — 1\ /n\ ( n \ 

~ (2n - ly- 1 + M*-1/W V(jfc + l)(n-fc + l)J 

n k fn\ ~^ d ~ 2 ^ f n \ 

~ (2n-l) d - 1 + ^ 2 n\k) \{k + l)(n - fc + 1) J 



Since fc/n < 1, n/(fe + l)(n — A; + 1) < 1/2, and (?) > f"), and there are fewer than 
n — 1 terms in the sum, the second term is at most 



n 



-(d-2) 



( n _l)-(d-3). 2 d-2. 2 -(d-l)/2 < 



1 



n 

as required. □ 



6. ASYMPTOTICS OF THE NUMBERS W n (d) 

Recall that w n (d) and u n (d) are the numbers of unrestricted, respectively indecom- 
posable, (i-dimensional 2-magic matrices of size n. Using the exponential principle, 
we can relate the sequence (w n (d)) n >o to the sequence (u n {d)) n >o for each fixed d, see 
( 16. ip below. This relationship combined with the fact that the sequence (u n (d)) n >o 
grows sufficiently rapidly for d > 3 (Proposition H] says that it grows very roughly like 
((2n)!) d_1 ) allows us to conclude that, for d > 3, w n (d) and u n (d) grow at the same 
rate. 

Theorem 5. For fixed d>3, we have 

w n {d) ~ 2- nd ((2n)\) d -\ asn^oo. 

Proof. Generalising the argument at the end of Section [3j we observe that d- dimensional 
2-magic matrices may be seen as a d-sort species in the sense of Joyal [T4] (see also [5J 
Def. 4 on p. 102]), with the row indices and the column indices forming the two set on 



which the functor defining the species operates. Hence, by the exponential principle for 
(i-sort species [TU Prop. 20] (see also P, Sec. 2.4]), we have 



J2 W ^Jri)d = eX P 

n>0 ^ \n>l 

If we now differentiate both sides of this equation with respect to z and subsequently 
multiply both sides by z, then we obtain 



(n\r \" — ' ' ' (n\r I " \ " — ' ' ' \n\) d 

a n \n>l v ' / \n>l 

\ ^ . z n 



v n>l ^ / \n>0 

Comparison of coefficients of z n on both sides then leads to the relation 

w n (d) = ttn(ri) + X, ~ ( , j u k (d)w n _ k (d). (6.1) 

As we said at the beginning of this section, our goal is to show that w n (d) grows 
asymptotically at the same rate as u n (d). Hence, putting w n (d) = u n (d) + x n (d), we 
have to show that x n (d) = o(u n (d)). We assume inductively that 

x m (d) < 2~ m ((2m- l)!!) d - 1 (m!) d - 1 

for all m between 2 and n — 1; the induction starts since we have X\(d) = X2{d) = 0. 
Now, using the inductive hypothesis with the recurrence relation (16. ID . we have 



((2n- l)!!)^ 1 ^!)^ 1 ~ ^n\kj VW/ 



< 



n ~ l fc /n\ - {d - 2) ( n n (d - 1)/2 



2 V^ 

n 



i n\kJ \(k + l)(n - fc + 1) 



< (2 1 / 2 n)-( d " 3 ), 

which establishes the result if d > 3. For d = 3, this inequality gives the inductive step 
(that is, that the left-hand side is at most 1); the fact that it is o(l) for large n is proved 
by an argument like that in the proof of Proposition |H □ 

Acknowledgment. The authors are grateful to a referee for some helpful information on 
contingency tables. 

References 

[1] M. Ahmed, J. A. De Loera and R. Hemmecke, Polyhedral cones of magic cubes and squares, in: 
Discrete and Computational Geometry. The Goodman-Pollack Festschrift (B. Aronov, S. Basu, 
J. Pach and M. Sharir, eds.), Springer- Verlag, Berlin, Algorithms Comb. 25, 2003, pp. 25-41. 

[2] M. Aigner, Combinatorial Theory, Springer- Verlag, Berlin, 1979. 

[3] H. Anand, V. C. Dumir, and H. Gupta, A combinatorial distribution problem, Duke Math. J. 33 
(1966), 757-769. 

[4] M. Beck, M. Cohen, J. Cuomo and P. Gribelyuk, The number of "magic" squares, cubes, and 

hypercubes, Amer. Math. Monthly 110 (2003), 707-717. 
[5] G. Birkhoff, Tres observaciones sobre el algebra lineal, Univ. Nac. Tucum'an Rev. Ser. A 5 (1946), 

147-150. 



A NOTE ON HIGHER-DIMENSIONAL MAGIC MATRICES 



9 



F. Bergeron, G. Labelle, and P. Leroux, Combinatorial Species and Tree-Like Structures, Cam- 
bridge University Press, Cambridge, 1998. 

P. J. Cameron, C. Krattenthaler and T. W. Miiller, Decomposable functors and the exponential 
principle, II, Seminaire Lotharingien Combin. 61A (2011), Article B6lAm, 38 pages. 
J. A. De Loera, F. Liu and R. Yoshida, A generating function for all semi-magic squares and the 
volume of the Birkhoff polytope, J. Algebraic Combin. 30 (2009), 113-139. 

J. A. De Loera and S. Onn, The complexity of three-way statistical tables, SIAM J. Comput. 33 
(2004), 819-836 (electronic). 

J. A. De Loera and S. Onn, All linear and integer programs are slim 3-way transportation pro- 
grams, SIAM J. Optim. 17 (2006), 806-821 (electronic). 

P. Diaconis and A. Gangolli, Rectangular arrays with fixed margins, in Discrete Probability and 
Algorithms (D. Aldous, P. Diaconis, J. Spencer, and J. M. Steele (eds.), IMA Vol. Math. Appl., 
Springer, New York, pp. 15-41. 

P. Flajolet and R. Sedgewick, Analytic Combinatorics, Cambridge University Press, Cambridge, 
2009. 

S. Hosten and S. Sullivant, Grdbner bases and polyhedral geometry of reducible and cyclic models, 
J. Combinatorial Theory Ser. A 100 (2002), 277-301. 

A. Joyal, Une theorie combinatoire des series formelles, Adv. in Math. 42 (1981), 1-82. 
P. A. MacMahon, Combinatory Analysis, vol. 2, Cambridge University Press, 1916; reprinted by 
Chelsea, New York, 1960. 

L. J. Slater, Generalized Hypergeometric Functions, Cambridge University Press, Cambridge, 1966. 
N. J. A. Sloane (editor), The On-Line Encyclopedia of Integer Sequences, available at 
http : //www. research. att . com: 80/~nj as/ sequences/. 

R. P. Stanley, Enumerative Combinatorics, vol. 1, Wadsworth & Brooks/Cole, Pacific Grove, 
California, 1986; reprinted by Cambridge University Press, Cambridge, 1998. 

B. Sturmfels, Grdbner Bases and Convex Polytopes, University Lecture Series 8, Amer. Math. 
Soc, Providence, RI, 1996. 

School of Mathematical Sciences, Queen Mary & Westfield College, University of 
London, Mile End Road, London El 4NS, United Kingdom. 

WWW: http : //www . maths . qmw . ac . uk/~p j c/ , http : //www . maths . qmw . ac . uk/~twm/ . 



Fakultat fur Mathematik, Universitat Wien, Nordbergstrasse 15, A-1090 Vienna, 
Austria. WWW: |http : //www .mat .univie . ac . at/-kratt|. 



